Paediatric leukaemia DNA methylation profiling using MBD enrichment and SOLiD sequencing on archival bone marrow smears

BACKGROUND: Acute Lymphoblastic Leukaemia (ALL) is the most common cancer in children. Over the past four decades, research has advanced the treatment of this cancer from a less than 60% chance of survival to over 85% today. The causal molecular mechanisms remain unclear. Here, we performed sequencing-based genomic DNA methylation profiling of eight paediatric ALL patients using archived bone marrow smear microscope slides. FINDINGS: SOLiD™ sequencing data was collected from Methyl-Binding Domain (MBD) enriched fractions of genomic DNA. The primary tumour and remission bone marrow sample was analysed from eight patients. Four patients relapsed and the relapsed tumour was analysed. Input and MBD-enriched DNA from each sample was sequenced, aligned to the hg19 reference genome and analysed for enrichment peaks using MACS (Model-based Analysis for ChIP-Seq) and HOMER (Hypergeometric Optimization of Motif EnRichment). In total, 3.67 gigabases (Gb) were sequenced, 2.74 Gb were aligned to the reference genome (average 74.66% alignment efficiency). This dataset enables the interrogation of differential DNA methylation associated with paediatric ALL. Preliminary results reveal concordant regions of enrichment indicative of a DNA methylation signature. CONCLUSION: Our dataset represents one of the first SOLiD™MBD-Seq studies performed on paediatric ALL and is the first to utilise archival bone marrow smears. Differential DNA methylation between cancer and equivalent disease-free tissue can be identified and correlated with existing and published genomic studies. Given the rarity of paediatric haematopoietic malignancies, relative to adult counterparts, our demonstration of the utility of archived bone marrow smear samples to high-throughput methylation sequencing approaches offers tremendous potential to explore the role of DNA methylation in the aetiology of cancer.

Genome sequencing: a grin without a cat

In some types of unicellular algae, the chloroplasts have their own nucleus — a legacy of the time when the chloroplast was a free-living cell. The sequence of the genome in one such nucleus is now revealed.

Eulerian superpath approach to correct sequencing error in shotgun assembly

This work describes an error correction method based on the Euler Superpath problem. Sequence data is mapped to an Euler Superpath dynamically by Merging Transformation. With restriction and guiding rules, data consistency is maintained and error paths are separated from correct data: Error edges are mapped to the correct ones and after substitution (of error edges with right paths), corresponding errors in the sequencing data are eliminated.

Modeling a just-in-time plant using sequencing algorithms

An Australian automotive component company plans to assemble and deliver seats to customer on just-in-time basis. The company management has decided to model operations of the seat plant to help them make decisions on capital investment and labour requirements. There are four different areas in seat assembly and delivery areas. Each area is modeled independently to optimise its operations. All four areas are then combined into one model called the plant model to model operations of seat plant from assembly to delivery. Discrete event simulation software is used to model the assembly operations of seat plant.

A pHMM-ANN based discriminative approach to promoter identification in prokaryote genomic contexts

The computational approach for identifying promoters on increasingly large genomic sequences has led to many false positives. The biological significance of promoter identification lies in the ability to locate true promoters with and without prior sequence contextual knowledge. Prior approaches to promoter modelling have involved artificial neural networks (ANNs) or hidden Markov models (HMMs), each producing adequate results on small scale identification tasks, i.e. narrow upstream regions. In this work, we present an architecture to support prokaryote promoter identification on large scale genomic sequences, i.e. not limited to narrow upstream regions. The significant contribution involved the hybrid formed via aggregation of the profile HMM with the ANN, via Viterbi scoring optimizations. The benefit obtained using this architecture includes the modelling ability of the profile HMM with the ability of the ANN to associate elements composing the promoter. We present the high effectiveness of the hybrid approach in comparison to profile HMMs and ANNs when used separately. The contribution of Viterbi optimizations is also highlighted for supporting the hybrid architecture in which gains in sensitivity (+0.3), specificity (+0.65) and precision (+0.54) are achieved over existing approaches.

CIDB: Chlamydia Interactive Database for cross-querying genomics, transcriptomics and proteomics data

Chlamydiae are important pathogens of humans, birds and a wide range of animals. They are a unique group of bacteria, characterized by their developmental cycle. Chlamydia has been difficult to study because of their obligate intracellular growth habit and lack of a genetic transformation system. However, the past 5 years has seen the full genome sequencing of seven strains of Chlamydia and a rapid expansion of genomic, transcriptomic (RT-PCR, microarray) and proteomic analysis of these pathogens. The Chlamydia Interactive Database (CIDB) described here is the first database of its type that holds genomic, RT-PCR, microarray and proteomics data sets that can be cross-queried by researchers for patterns in the data. Combining the data of many research groups into a single database and cross-querying from different perspectives should enhance our understanding of the complex cell biology of these pathogens. The database is available at: http://www3.it.deakin.edu.au:8080/CIDB/.

Genomic analyses and cloning of novel chicken natriuretic peptide genes reveal new insights into natriuretic peptide evolution

The natriuretic peptide (NP) family consists of multiple subtypes in teleosts, including atrial, B-type, ventricular, and C-type NPs (ANP, BNP, VNP, CNP-1–4, respectively), but only ANP, BNP, CNP-3, and CNP-4 have been identified in tetrapods. As part of understanding the molecular evolution of NPs in the tetrapod lineage, we identified NP genes in the chicken genome. Previously, only BNP and CNP-3 have been identified in birds, but we characterized two new chicken NP genes by cDNA cloning, synteny and phylogenetic analyses. One gene is an orthologue of CNP-1, which has only ever been reported in teleostei and bichir. The second gene could not be assigned to a particular NP subtype because of high sequence divergence and was named renal NP (RNP) due to its predominant expression in the kidney. CNP-1 mRNA was only detected in brain, while CNP-3 mRNA was expressed in kidney, heart, and brain. In the developing embryo, BNP and RNP transcripts were most abundant 24 h post-fertilization, while CNP mRNA increased in a stage-dependant manner. Synthetic chicken RNP stimulated an increase in cGMP production above basal level in chicken kidney membrane preparations and caused a potent dose-dependant vasodilation of pre-constricted dorsal aortic rings. From conserved chromosomal synteny, we propose that the CNP-4 and ANP genes have been lost in chicken, and that RNP may have evolved from a VNP-like gene. Furthermore, we have demonstrated for the first time that CNP-1 is retained in the tetrapod lineage.

A genomic approach to understanding the cause and effect of annual ryegrass toxicity

Annual Ryegrass Toxicity (ARGT) is a potentially lethal disease affecting livestock grazing on pastures or consuming fodder that include annual ryegrass (Lolium rigidum) contaminated with corynetoxins. The corynetoxins (CTs), among the most lethal toxins produced in nature, are produced by the bacterium Rathayibacter toxicus that uses a nematode vector to attach to and infect the seedheads of L.rigidum. There is little known of the factors that control toxin production. Several studies have speculated that a bacteriophage specific to R.toxicus may be implicated in CT production. We have developed a PCR-based assay to test for both bacterium and phage in ryegrass material and results indicate that there is a correlation between phage and bacterial presence in all toxic ryegrass samples tested so far. This PCR-based technique may ultimately allow for a rapid, high-throughput screening assay to identify potentially toxic pastures and feed in the field. Currently, ~80% of the 45 Kb genome has been sequenced an investigation to further elucidate its potential role in toxin production.Furthermore, specific alterations in gene expression as a result of exposure to CTs or the closely related tunicamycins (TMs), which are commercially available and considered biologically indistinguishable from CTs, will be evaluated for use as biomarkers of exposure. The effects of both toxins will be analysed in vitro using a rat hepatocyte cell line and screened on a low-density DNA micro array “CT-Chip” that contains <100 selected rat hepatic genes. The results are expected to further define the bioequivalence of CTs and TMs and to identify levels of exposure that are related to specific toxic effects or have no adverse effect on livestock.

Virulence determinants of infectious bursal disease virus

The very virulent (vv) pathotype of infectious bursal disease virus (IBDV) has spread rapidly throughout Europe, Asia, and the Middle East. Although Australia is currently unaffected, there remains the potential for incursion of an exotic isolate. The aim of this study was to identify putative virulence determinants of IBDV to facilitate the development of improved diagnostic assays for detection and characterisation of vvIBDV isolates. Sequencing of Indonesian vvIBDV Tasik94 revealed a unique substitution [ A�¨S222] in the hypervariable region (HVR) of viral protein (VP) VP2, which did not appear to impinge on virulence or antigenicity. Phylogenetic analyses indicated that Tasik94 was closely related to Asian and European vvIBDV strains. Extensive alignment of deduced protein sequences across the HVR of VP2 identified residuesI242 I256 and I294 as putative markers of the vv phcnotype. Comparison of the pathology induced by mildly-virulent Australian IBDV 002/73 and Indonesian vvIBDV Tasik94, revealed that histological lesions in the spleen, thymus and bone marrow were restricted to Tasik94-infected birds, suggesting the enhanced pathogenicity of vvIBDV might be attributed to replication in non-bursal lymphoid organs. The biological significance of the VP2 HVR in virulence was assessed using recombinant viruses generated by reverse genetics. Both genomic segments of Australian IBDV 002/73, and recombinant segment A constructs in which the HVR of 002/73 was replaced with the corresponding region of either tissue culture-adapted virus or vvIBDV (Tasik94), were cloned behind T7 RNA polymerase promoter sequences. In vitro transcription/translation of each construct resulted in expression of viral proteins. Co-transfection of synthetic RNA transcripts initiated replication of both tissue culture-adapted parental and recombinant viruses, however attempts to rescue non-adapted viruses in specific-pathogen-free (SPF) chickens were unsuccessful. Nucleotide sequence variation in the HVR of VP2 was exploited for the development of a new diagnostic assay to rapidly detect exotic IBDV isolates, including vvIBDV, using reverse transcription polymerase chain reaction (RT-PCR) amplification and Bmrl restriction enzyme digestion. The assay was capable of differentiating between endemic and exotic IBDV in 96% of 105 isolates sequenced to date.

Familial and genomic analyses of postural changes in systolic and diastolic blood pressure

The physiological adaptation to the erect posture involves integrated neural and cardiovascular responses that might be determined by genetic factors. We examined the familial- and individual-specific components of variance for postural changes in systolic and diastolic blood pressure in 767 volunteer nuclear adult families from the Victorian Family Heart Study. In 274 adult sibling pairs, we made a genome-wide scan using 400 markers for quantitative trait loci linked with the postural changes in systolic and diastolic pressures. Overall, systolic pressure did not change on standing, but there was considerable variation in this phenotype (SD=8.1 mm Hg). Familial analyses revealed that 25% of the variance of change in systolic pressure was attributable to genetic factors. In contrast, diastolic pressure increased by 6.3 mm Hg (SD=7.0 mm Hg) on standing and there was no evidence of contributory genetic factors. Multipoint quantitative genome linkage mapping suggested evidence (Z=3.2) of linkage of the postural change in systolic pressure to chromosome 12 but found no genome-wide evidence of linkage for the change in diastolic pressure. These findings suggest that genetic factors determine whether systolic pressure decreases or increases when one stands, possibly as the result of unidentified alleles on chromosome 12. The genetics of postural changes in systolic blood pressure might reflect the general buffering function of the baroreflex; thereby, the predisposition to sudden decreases or increases in systolic pressure might cause postural hypotension or vessel wall disruption, respectively.